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I. INTRODUCTION 



Quantum-classical analogies have been exploited on 
many occasions to mimic at a macroscopic level many 
quantum phenomena which are currently inaccessible in 
microscopic quantum systems In particular, in the 
past decade engineered photonic structures have pro- 
vided a useful laboratory tool to investigate and visu- 
alize with classical optics the dynamical aspects em- 
bodied in a wide variety of coherent quantum phenom- 
ena encountered in atomic, molecular, condensed-matter 
and matter- wave physics 0. Among others, we men- 
tion the optical analogues of electronic Bloch oscillations 
3 , H[ and Zener tunneling Q , dynamic localization 
7|, coherent enhancement and destruction of tunneling 
8j|, adiabatic stabilization of atoms in ultrastrong laser 
fields 3, Anderson localization [Hj, quantum Zeno ef- 
fect [Tl|, coherent population transfer and coher- 
ent vibrational dynamics [Hj]. Most of the above men- 
tioned studies are based on the formal analogy between 
paraxial wave equation of optics in dielectric media and 
the single-particle nonrelativistic Schrodinger equation, 
thus providing a test bed for classical analogue stud- 
ies of nonrelativistic quantum phenomena. Recently, a 
great attention has been devoted toward the investiga- 
tion of experimentally-accessible and controllable classi- 
cal or quantum systems that simulate certain fundamen- 
tal phenomena rooted in the relativistic Dirac equation. 
Among others, cold trapped atoms, ions and graphene 
have proven to provide accessible systems to simulate 
relativistic physics in the lab, and a vast literature on 
this subject has appeared in the past few years (see, for 
instance, fl4l - flr| | and references therein). In particular, 
low-energy excitations of nonrelativistic two-dimensional 
electrons in graphene obey the Dirac- Weyl equation and 
behave like massless relativistic particles. This has lead 
to the predictions in condensed-matter physics of phe- 
nomena analogous to Zitterbewegung [I?} and Klein tun- 
neling 18] of relativistic massive or massless particles, 
with the first experimental evidences of Klein tunneling 
in graphene [l9[ and in carbon nanotube [2(| systems. 
Photonic analogues of Dirac-type equations have been 



also theoretically proposed for light propagation in cer- 
tain triangular or honeycomb photonic crystals, which 
mimic conical singularity of energy bands of graphene 
[2ll [22j . as well as in metamaterials [23[ and in opti- 
cal superlattices [24j. This has leads to the proposals of 
photonic analogues of relativistic phenomena like Zitter- 
bewegung [22l - [24j and Klein tunneling [2^, [25[ . 
Electron-positron pair production due to the instability 
of the quantum electrodynamics (QED) vacuum in an 
external electric field (a phenomenon generally referred 
to as vacuum decay) is another remarkable prediction of 
Dirac theory and regarded as one of the most intrigu- 
ing non-linear phenomena in QED (see, for instance, 
[2a . [27(| ) . In intuitive terms and in the framework of 
the one-particle Dirac theory, the pair production process 
can be simply viewed as the transition of an electron of 
the Dirac sea occupying a negative-energy state into a fi- 
nal positive-energy state, leaving a vacancy (positron) in 
the negative-energy states. There are basically two dis- 
tinct transition mechanisms: the Schwinger mechanism 
in presence of an ultrastrong static electric field, and dy- 
namic pair creation in presence of time- varying electric 
fields. The Schwinger mechanism J28] can be understood 
as a tunneling process through a classically forbidden re- 
gion, bearing a close connection to Klein tunneling. Dy- 
namic pair creation, originally proposed by Brezin and 
Itzykson [29| for oscillating spatially-homogeneous fields 
and subse quently investigated by several authors (see, for 
instance, |30l - [32| and references therein), has attracted 
recently great attention since the advent of ultrastrong 
laser system facilities, which pave the way toward the 
realization of purely laser-induced pair production [33| . 
Electric fields oscillating in time only can be achieved at 
the antinodes of two oppositely propagating laser beams, 
and can lead to such intriguing phenomena as Rabi os- 
cillations of the Dirac sea (see, for instance, [27|,|3l|). In 
the framework of the one-particle Dirac theory of vac- 
uum decay [3lj |. a simple picture of pair production is 
represented by the time evolution of an initially negative- 
energy Gaussian wave packet, representing an electron in 
the Dirac sea, under the influence of an oscillating elec- 
tric field [HJ . When the e + e~ pair is produced, a droplet 
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is separated from the wave packet and moves opposite to 
the initial one [32[ . The droplet is a positive-energy state 
and represents the created electron. Such a dynamical 
scenario of pair production, observed in numerical simu- 
lations [33|, is unlikely to be accessible in a foreseen ex- 
periment using two counterpropagating and ultraintense 
laser beams, therefore it may be of some interest to find 
either a classical or a quantum simulator capable of visu- 
alizing in the lab the wave packet dynamics of e + e~ pair 
creation. It is the aim of this work to propose a classic 
wave optics analogue of the QED e + e~ pair production 
in oscillating fields, based on monochromatic light prop- 
agation in curved waveguide optical superlattices [g, [34| , 
in which the temporal evolution of Dirac wave packets 
and e + e _ pair production is simply visualized as spatial 
beam break up along the curved photonic structure. As 
our analysis is focused onto a classical wave optics ana- 
logue of QED vacuum decay, a similar dynamical scenario 
could occur for ultracold atoms in accelerated double- 
periodic optical lattices (see, for instance, [Hj]), which 
might thus provide a quantum simulator for relativistic 
QED decay complementary to our classical analogue. 



II. BASIC MODEL AND QUANTUM-OPTICAL 
ANALOGY 

A. The photonic structure 

The photonic structure considered in this work con- 
sists of a one-dimensional binary supcrlatticc of weakly- 
guiding dielectric optical waveguides with a periodically- 
curved optical axis. This optical system was previously 
introduced in another physical context for multiband 
diffraction and refraction control of light waves in optical 
lattices [34j. It should be noted that a similar photonic 
structure, but with a circularly-curved optical axis, has 
been recently realized to visualize in the lab the optical 
analogue of coherent Bloch-Zener oscillations of nonrel- 
ativistic electrons in crystalline potentials @, whereas 
the same structure but with a straight optical axis has 
been proposed to realize a photonic analogue of Zitterbe- 
wegung and Klein tunneling of relativistic electrons [24| . 
The photonic superlattice consists of two interleaved and 
equally spaced lattices A and B, as shown in Fig. 1(a), 
with spacing a between adjacent waveguides. The nor- 
malized waveguide index profile is assumed to be the 
same for the waveguides of the two lattices A and B, 
with an alternation of peak refractive index changes dn\ 
and dn2 for the two lattices [6|. The waveguide axis is 
periodically bent along the propagation direction z [see 
Fig. 1(b)], with an axis bending profile Xq(z) of period A; 
typically the condition i(0) = is assumed to ensure that 
the waveguide axis is orthogonal to the exciting input 
plane z — 0, where the dot denotes the derivative with 
respect to z. Light propagation in the curved superlat- 
tice is at best captured in the 'waveguide' reference frame 
after a Kramers-Henneberger transformation, where the 



array appears to be straight (see, for instance, Q)- In 
this reference frame, the waveguide curvature along the 
propagation axis appears as a fictitious refractive index 
gradient in the transverse plane with a z-varying slope 
proportional to the local waveguide axis curvature 
After a gauge transformation for the electric field enve- 
lope, propagation of monochromatic waves at wavelength 
(in vacuum) A is described by the following scalar wave 
equation (see, for instance, pi l34jl 

dE A d 2 E n s — n(x) 2Tm s x (z) „ 

ij- = — ~o 2" + 2lT I E + -T^xE 

oz 47m s ox z A A 

where n s is the bulk refractive index, nix) is the refrac- 
tive index profile of the superlattice, and the double dot 
stands for the second derivative with respect to z. In 
the tight-binding approximation, light transport in the 
curved photonic lattice can be described by means of 
coupled-mode equations for the amplitudes q of light 
modes trapped in the various waveguides [33 ] 

* J = -a[c l+ i + q_i) + (-l) l S Cl + F(z)lci, (2) 

where 26 and a are the propagation constant mismatch 
and the coupling rate between two adjacent waveguides of 
lattices A and B, and F(z) = 2irn s axo(z) / 'A. The tight- 
binding model (2) is accurate to describe beam dynamics 
of Eq.(l) provided that the first two minibands of the 
array are involved in the dynamics (34|. A similar tight- 
binding model also applies to the dynamics of cold atoms 
in accelerated double-periodic optical lattices (35|. F° r 
a straight lattice, i.e. for F(z) — 0, the tight-binding 
minibands of the superlattice are readily calculated by 
making the plane- wave Ansatz Ci(q) ~ exp(iqla — icoz) in 
Eq.(2) and read [see Fig. 1(c)] 

ui±(q) =±y/6 2 + 4ct 2 cos 2 (ga) (3) 
where q is the wave number. 

B. Quantum-optical analogy 

In Ref . [24j . it was recently shown that, near the Bril- 
louin zone edge q = ±7r/(2a), the dispersion relations 
(3) of the straight superlattice approximate the positive 
(electron) and negative (positron) energy curves of a one- 
dimensional massive Dirac electron in absence of exter- 
nal fields, and that beam propagation at incidence angles 
near the Bragg angle mimics the temporal dynamics of a 
one-dimensional free Dirac electron, the two components 
of the spinor wave function corresponding to the occu- 
pation amplitudes in the two sublattices A and B of the 
superlattice. To study the photonic analogue of vacuum 
decay and electron-positron production arising from an 
oscillating field, we consider here light propagation in a 
periodically-curved binary array within the tight-binding 
model (2). As it will be shown below, the undulation of 
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FIG. 1: (color online) (a) and (b) Schematic of a binary su- 
perlattice, made of two interleaved lattices A and B of waveg- 
uides with high (dni) and low (dri2) refractive index changes, 
equally spaced by a [Fig. 1(a)] and with a periodically-curved 
axis [Fig. 1(b)]. The array is excited by a broad beam tilted at 
the angle 8. (c) Dispersion curves (minibands) of the tight- 
binding binary superlattice [Eq.(3)] for the straight binary 
array (solid curves), and corresponding electron and positron 
dispersion curves of the Dirac equation (7) (in absence of the 
ac field) near the Brillouin zone edge q = n/(2a) (dotted 



the waveguide axis along the propagation direction mim- 
ics an external ac electric field in the one-dimensional 
Dirac equation. To explain the quantum-optical analogy, 
let us set in Eqs.(2) q(z) = ai(z) exp[— i<J>(z)Z], where 



/o A 
Coupled-mode equations than take the equivalent form 
.dai 



dz 



-a exp[— i$(z)]a; +1 — a exp[i<I>(z)]ai_i 



+ i-iysai. 



(5) 



Note that, for the special case of a sinusoidal axis bend- 
ing profile x (z) = —Acos(2ttz/ A) of amplitude A and 
period A, one has $>(z) = $0 sin(27rz/A) with 



$n = 



4n 2 n,aA 



AA 



(6) 



Let us then introduce the following two assumptions: (i) 
the amplitude of axis bending profile is small enough 
that 1 3? I <C 7r/2 over the entire oscillation cycle; (ii) the 
array is excited by a broad beam tilted at the Bragg 
angle 6b = A/(4n s a). The first assumption enables 
to set exp(±i$) ~ 1 ± i<f> in Eqs.(5); the second as- 
sumption implies that the modes in adjacent waveguides 
are excited with a nearly equal amplitude but with a 
phase difference of 7r/2 (24|. Therefore, after setting 
an{z) = (-l)Vi(Z, z), a 2 i-i = z), the am- 

plitudes ipi and ^2 vary slowly with I, and one can thus 
write ipi^Q ±l,z) — ipi,2{l,z) ± (dtpi^/dl) and con- 
sider I = £ = xj (2a) as a continuous variable rather 



than as an integer index. Under such assumptions, from 
Eqs.(5) it readily follows that the two-component spinor 
ip(£,z) = (ipi, "4>2) T satisfies the one-dimensional Dirac 
equation in presence of an external ac field, namely 



dz 



dtp 

—iaa— — 2(T<f>(z)a-0 + 5f3ip, 



where 



1 

1 



1 

-1 



(7) 



(8) 



are the a x and <j z Pauli matrices, respectively. Note that, 
after the formal change 



(9) 
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Eq.(7) corresponds to the one-dimensional Dirac equa- 
tion for an electron of mass m and charge e in presence of 
a spatially-homogeneous and time- varying vectorial po- 
tential A = (A x ,0,0), which describes the interaction 
of the electron with an external oscillating electric field 
E x (t) = — (dA x / dt) in the dipole approximation (see, 
for instance, [36|]). Note also that, in our optical ana- 
logue, the temporal evolution of the spinor wave func- 
tion ip for the Dirac electron is mapped into the spa- 
tial evolution of the field amplitudes ip\ and ip2 along 
the z-axis of the array, and that the two components 
"01! ^2 of the spinor wave function correspond to the 
occupation amplitudes in the two sublattices A and B 
composing the superlattice. In absence of the external 
field, the energy-momentum dispersion relation huj(k) of 
the Dirac equation (7), obtained by making the Ansatz 
ip ~ exp(«£;£ — iojz), is composed by the two branches 
w±(k) — ±e(fc), corresponding to positive and negative 
energy states of the relativistic free electron, where 



E(fc) = VS 2 +<7 2 k 2 . 



(10) 



Note that such two branch curves are readily obtained 
from Eq.(3) after setting 



q = 



2 a 



k 

2a. 



(11) 



and assuming small values of k, which corresponds to the 
assumption of slowly-varying envelopes tpi^l, z). There- 
fore, the positive and negative energy states of the free 
Dirac electron are mapped into the dispersion curves of 
the two minibands of the binary array near the boundary 
of the Brillouin zone at q = 7r/(2a), as shown in Fig. 1(c). 
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III. OPTICAL ANALOGUE OF PAIR 
PRODUCTION IN AN OSCILLATING FIELD 



The phenomenon of vacuum decay by oscillating fields, 
originally predicted in Ref.[2!|, refers to the electron- 
positron pair production due to the instability of the 
QED vacuum in a superstrong field of two counterprop- 
agating laser beams [3l|, [HI ■ Vacuum decay should be 
generally investigated in the framework of a quantum 
field theory of the Dirac equation, however a simpler 
and more intuitive approach is possible within the one- 
particle Dirac equation (7) , where pair production corre- 
sponds to a field-induced transition of an electron in the 
Dirac sea of neg ative-energy states to a final positive- 
energy state |3ll. |32|. In this approach, a simple picture 
of pair production as a wave packet break-up process 
has been suggested [32|. Here an electron in the nega- 
tive sea is represented by a Gaussian wave packet formed 
by a superposition of negative-energy states. When the 
e + e~ pair is produced under the action of the oscillating 
field, 



a droplet is separated from the wave packet and 
moves opposite to the initial one (see, for instance, Fig.l 
of Ref.[32]). The droplet is a positive-energy state and 
represents the created electron. It has been shown that 
such a simplified model of pair production, besides to be 
rather intuitive, also provides a correct way to calculate 
the rate of particle creation [32] . Our classical wave op- 
tics analogue, grounded on the formal analogy between 
the one-particle Dirac equation and the spinor-like wave 
equation (7) via the variable transformation (9), thus en- 
ables to visualize in space the analogue of pair produc- 
tion as a break up of an initial Gaussian wave packet, 
composed by a superposition of Bloch modes of the low- 
est lattice miniband and representing an electron in the 
Dirac sea, under the influence of an oscillating electric 
field. 



A. Two-level model 

For a spatially-homogeneous and time-dependent field, 
it is known that momentum conservation reduces the 
problem of pair creation to a two-level system consisting 
of a negative and a positive energy state coupled by the 
external field. Within the two-level model, pair produc- 
tion generally occurs as a multiphoton resonance process 
enforced by energy conservation, with interesting effects 
such as Rabi oscillations of the quantum vacuum (see, for 
instance, [3l|, HH). The two- level description of pair pro- 
duction in the dipole approximation for the one-particle 
Dirac equation can be found, for instance, in Refs.[3lj|; 
for the sake of completeness, it is briefly reviewed in the 
Appendix for the case of the one-dimensional Dirac equa- 
tion. Owing to the quantum-optical analogy established 
in Sec. II. B, the same dynamical scenario is thus expected 
for light transport in a periodically-curved tight-binding 
binary array. To this aim, let us look for a solution to 



coupled-mode equations (5) of the form 

si(z) 



ai{z) 



s 2 (z) 



exp(iqla) 



(12) 



where the wave number q is related to the electron mo- 
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FIG. 2: (color online) Transition probability P — \r+(z)\ 2 
versus propagation distance z, normalized to the ac oscillation 
period A, in a tight-binding binary lattice with a sinusoidally- 
curved axis for parameter values a = 2, 8 = 1.817, qa — tv/4, 
$0 = 0.4 and for two values of modulation period A: (a) A = 
2.8556, and (b) A = 2.7196. In (a) the multiphoton resonance 
condition (22) is satisfied for n = 3, and Rabi flopping is 
clearly observed. In (b) the modulation period is slightly 
decreased by 5% from the n = 3 resonance value. 

mentum k for the Dirac equation (7) by Eq.(ll), and 
where the upper (lower) row applies to an even (odd) 
value of index I. Substitution of Eq.(12) into Eq.(5) yields 
the following coupled equations 



dz 



T(z) 



(13) 



where 



T(z) = 



-2<rcos[qa - <&(z)] 



-2a cos[ga — $>(z) 

-S 



(14) 

For a straight array, i.e. in absence of the external 
field $ = 0, the (normalized) eigenvectors of the z- 
independent matrix T, 



2a cos(qa) 



u±(q)-8 

y/2\v±(q)\\u±(q)-6\ 



(15) 



corresponding to the eigenvalues u±(q) defined by Eq.(3), 
belong to the positive (electron) and negative (positron) 
energy branches. These states are basically the Bloch 
modes of the superlattice, belonging to the first two low- 
est bands, with wave number q. If the initial state at 
z = is prepared in the Bloch mode belonging to the 
negative energy branch (corresponding to an electron oc- 
cupying a negative energy state of the Dirac sea in the 
quantum mechanical context), i.e. if (si, S2) T = (si, ^2)^ 
at z = 0, the application of the oscillating external field 
$(z) can induce a transition to the Bloch mode (si, S2) + 
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on the positive energy branch, i.e. the optical analogue of 
an electron-positron pair production is realized. It should 
be noted that direct (and even indirect) interband tran- 
sitions and Rabi oscillations induced by a suitable lon- 
gitudinal refractive index modulation in singly-periodic 
waveguide arrays have been recently predicted and ex- 
perimentally observed in Refs.[37|. however in those pre- 
vious works the transitions are not described in terms of 
a Dirac-type equation with an oscillating field, and the 
quantum-optical analogy is thus not evident [38| ■ 
To study the field-induced transition process, it is worth 
projecting the state (si,S2) onto the Bloch mode basis 
(15), i.e. let us set 



-(*) 



r+(z) 



(16) 



where r—(z) and r + (z) are the occupation amplitudes 
of the negative and positive electron branches, respec- 
tively. Substitution Eq.(16) into Eq.(13) yields the fol- 
lowing coupled equations for the occupation amplitudes 



dz 



r_i_ 



= Z(z) 



(17) 



where the z-varying elements of the 2x2 matrix Z(z) 
are given by 

S 2 + 4(T 2 cos(ga) cos[ga - $(z)} 

-£■11 — —^22 — — \ L °) 



Z\2 — Zi\ 



2a5 {cos(qa) — cos[qa — $>(z)]} 



(19) 



Equations (17) are generally encountered in problems of 
driven two-level systems, and therefore phenomena like 
population Rabi flopping, multiphoton resonances, etc. 
are expected to occur in our system as well. Note that 
the conservation of population implies that |r_(z)| 2 + 
|r + (z)| 2 = 1. For the initial condition r_(0) = 1 and as- 
suming that the external ac field is switched on at z = 
and off after a propagation distance z = Lq (which is as- 
sumed to be an integer multiple of the modulation cycle 
A), i.e. &(z) =0 for z > Lo, the transition probability 
is thus given by P = |7-+(Lo)| 2 - The connection between 
the two- level model (17), derived from the tight-binding 
lattice model (5), and the two-level model of pair produc- 
tion in the framework of the one-dimensional Dirac equa- 
tion (7) is discussed in the Appendix. The latter model is 
retrieved from the former one by assuming a wave num- 
ber q near the edge of the Brilluoin zone [Eq.(ll) with 
|fc| <C 7r/2] and a small value of |3>(z)|; in this case, the 
expressions of matrix coefficients take the following sim- 
plified form [see also Eqs.(A9) and (A10) given in the 
Appendix] 



Zu = ~Z 22 = e(fe) - 



2a 2 k<^(z) 



2er<5$(z) 
^(fcT 



(20) 
(21) 



where e(k) is defined by Eq.(10). It is worth noticing 
that, under such an assumption, Eqs.(17) are analogous 
to the two-level equations describing population dynam- 
ics in a two-level dipolar molecule under an external ac 
field (see, for instance, |39(). 

In the context of field-induced pair production, the pho- 
ton energy 2ttH/ A of the ac field is generally much smaller 
that the energy separation 2mc 2 between negative and 
positive energy states, and the transition of one electron 
out of the Dirac sea into a positive-energy state gener- 
ally occurs via a multiphoton resonant process (see, for 
instance, (3lL l32l|). The condition of n-photon resonance, 
obtained in the low ac frequency limit by a WKB analysis 
of Eqs. (17-19), reads explicitly (see, e.g., [H]) 



2tt 



(22) 



where the quasi-energy £ (q) is calculated as 



£(q) = j J dz^Zf^ + ZUz) (23) 

= — / dy^/S 2 + Act 2 cos 2 [qa - $(Az//27r)]. 
27r Jo 

As an example, Fig. 2(a) depicts the behavior of the 
transition probability P = |r + (z)| 2 versus z, showing 
multiphoton Rabi oscillations, as obtained by numeri- 
cal simulations of Eqs. (17-19) for a sinusoidal ac field 
<&(z) = $o sin(27rz/A) for parameter values a = 2, 
S = 1.817, qa = 7r/4, $o = 0.4 and for a modula- 
tion period A = 2.8556, which satisfies the multiphoton 
resonance condition (22) with n = 3. For comparison, 
Fig. 2(b) shows the behavior of the transition probabil- 
ity for the same parameter values, but for a modulation 
period A = 2.7196 which is 5% smaller than the reso- 
nant value of Fig. 2(a). A comparison of Figs. 2(a) and 
(b) clearly indicates that, as it is well known, multipho- 
ton pair production is a resonant process. 
In our photonic system, the constraints of low modula- 
tion frequency and low field amplitudes typical of laser- 
driven QED vacuum can be overcome. For typical geo- 
metrical settings which apply to binary waveguide arrays 
and for reasonable propagation lengths achievable with 
current samples, the regimes of fast modulation frequen- 
cies (i.e., short bending periods A) and strong amplitudes 
are indeed more accessible than the low-frequency mul- 
tiphoton resonance regime; on the other hand, such a 
latter regime mi ght be accessible for cold atoms in opti- 
cal superlattices [35[. For instance, the values of a and 5 
used in Fig. 2 correspond to the coupling rate and propa- 
gation constant mismatch, in units of cm -1 , of a typical 
binary waveguide array for parameters values discussed 
in the next subsection. For such an array, the modula- 
tion period of axis bending in Fig. 2(a) that achieves the 
n = 3 photon resonance condition is A = 2.8556 cm, 
and the multiphoton Rabi oscillation period shown in 
Fig. 2(a) would thus correspond to about 80 cm, a length 
which is not accessible in an experiment. Conversely, the 
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regimes of fast modulation frequency and strong bending 
amplitudes are easily accessible. In such regimes, efficient 
transition can occur even for a single-cycle of the ac field 
[$(z) = $ sin(27rz/A) for < z < A, $(z) = for 
z < and z > A] , corresponding to the use of ultrastrong 
single-cycle counter-propagating laser pulses in the QED 
context. As an example, Fig. 3(a) shows the behavior of 
the transition probability P, after the interaction with a 
single-cycle pulse, as obtained by numerical simulations 
of Eqs.(17-19) for a = 2, S = 1.817, qa = tt/4 (as in 
Fig. 2) but with a much shorter period A = 0.6676, and 
for increasing values of the the field amplitude $o- Ex- 
amples of the detailed behavior of P(z) = |r + (z)| 2 along 
the oscillation cycle, corresponding to $o = 4 and $o = 6 
[points A and B in Fig. 3(a)] , are depicted in Figs. 3(b) and 
(c), respectively. 




amplitude <t* z /\ z /a 



FIG. 3: (color online) Transition probability P in a tight- 
binding binary lattice with a single-cycle sinusoidally-curved 
axis for parameter values a = 2, 8 — 1.817, qa = 7r/4, 
A = 0.6676, and for increasing values of amplitude 4>o. Points 
A, B and C in the figure correspond to the values 3>o — 4, 
$o = 6 and $0 = 0, respectively, (b) and (c) show the de- 
tailed behavior of |r+(z)| 2 versus propagation distance z in 
the single-cycle modulation period corresponding to points A 
[Fig.3(b)] and B [Fig.3(c)] of Fig.3(a). 



B. Wave packet dynamics 

The two-level description of pair production in a 
spatially-homogeneous and oscillating field described in 
the previous subsection assumes states with definite mo- 
mentum p in the negative and positive energy branches, 
and thus fully delocalized in space. A better visualization 
of the pair creation process is attained by considering, as 
an initial state, a free wave packet in the negative-energy 
continuum (for instance Gaussian-shaped), representing 
an electron in the Dirac sea [32]. After the external ac 
field is switched on for some interval (for instance for 
one oscillation cycle) and then switched off again, tran- 
sitions into the positive-energy continuum is visualized 
as a break up of the initial wave packet into two wave 
packets, which propagates with different group velocities 
and thus separate in space after some time [32[. These 
two wave packets represent the amplitude probabilities 
for the electron to be excited in the positive-energy con- 
tinuum or to remain in the Dirac sea. In our photonic 



analogue, wave packet break up is simply explained be- 
cause of the different refraction angles of wave packets 
belonging to the two minibands of the superlattices: the 
creation of a wave packet in the upper lattice miniband, 
induced by the longitudinal bending of the waveguide 
axis, is simply observed as a deviation of the propaga- 
tion direction from that of the initial wave packet. 
We have checked such a scenario by direct numerical sim- 
ulations of the paraxial field equation (1) using a stan- 
dard pseudospectral split-step method. Parameter values 
and refractive index profiles used in the simulations are 
compatible with binary wavegude arrays realized in fused 
silica by femtosecond laser writing and excited in the visi- 
ble at A = 633 nm @. Figure 4(a) shows the index profile 
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FIG. 4: (color online) (a) Refractive index profile of the bi- 
nary superlattice used in numerical simulations, and (b) cor- 
responding band diagram. Parameter values are given in the 
text. 

n(x) of the superlattice used in the simulations, corre- 
sponding to a waveguide spacing a = 10 /im, refractive 
index changes dn\ — 0.002 and drii = 0.00196 of adjacent 
waveguides, and a substrate refractive index n s = 1.42. 
Figure 4(b) shows the dispersion curves uj n (q) of a few 
low-order bands of the straight superlattice (band dia- 
gram). uj„(q) are computed by a standard plane- wave 
expansion method by looking for a solution to Eq.(l), 
with F = 0, of Bloch-Floquet type, i.e. of the form 
E(x, z) — u n {x, q) exp(igx) exp[—iui n (q)z], where n is the 
band order, q varies in the first Brillouin zone — n/ (2a) < 
q < 7r/ (2a), and u n (x, q) in the periodic part of the Bloch 
mode [u n {x + 2a, q) = u n {x, q)\. Note that the two lowest 
bands n = 1 and n = 2 in Fig. 4(b) correspond to the two 
minibands o;_ and uj + , respectively, of the tight-binding 
model (2) depicted in Fig. 1(c), i.e. to the negative- and 
positive-energy branches of the Dirac equation (7) in ab- 
sence of the external field. For the binary array of Fig. 4, 
the coupling rate a and propagation constant mismatch 
6 entering in the tight-binding model (2) can be simply 
estimated by fitting the two lowest bands of Fig. 4(b) us- 
ing Eq.(3) with a and S as fitting parameters, yielding 
a ~ 2 cm -1 and 5 ~ 1.817 cm -1 . To visualize the op- 
tical analogue of pair production, the bending profile of 
the waveguide axis is designed to simulate the action of a 
single-cycle ultrastrong field in the QED context, namely 
xq(z) — —Acos(2irz/ A) for < z < A, and Xq(z) = —A 
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for z > A. The array is excited at the input plane by a 
broad Gaussian beam of spot size wq, tilted at the an- 
gle 9, i.e. Eq.(l) is integrated with the initial condition 
E(x,0) — exp[— (x/wq) 2 ) exp(2irin s x9 / A) . The tilt angle 
has been chosen half of the Bragg angle, i.e. 9 = 9b /2, 
where 9b — 0.64°. At such a relatively small excita- 
tion angle, the lowest band of the array is mostly excited 
(see, for instance, [34[), and therefore the Gaussian wave 
packet mimics an initial electron wave packet of the Dirac 
sea, i.e. in the negative-energy spectrum. For 9 = 9 B /2, 
the condition qa — 7r/4 of Fig. 3 is also satisfied. The 
modulation period has been set equal to A ~ 0.67 cm to 
reproduce the conditions of Fig. 3. Figure 5 shows a few 
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FIG. 5: (color online) Beam propagation in a 6-cm-long one- 
dimensional binary waveguide array (snapshot of \E(x, z)\ 2 ), 
comprising a first section with a single-cycle sinusoidally 
curved axis of period A = 6.67 and amplitude A, and a sec- 
ond section with a straight axis: (a) A — (straight array), 
(b) A = 30 fim, and (c) A = 45 fim The array is excited by 
a Gaussian beam tilted at half of the Bragg angle (parameter 
values are given in the text). 

examples of beam propagation along a 6-cm-long array 
for an initial beam spot size wq — 80 /im and for increas- 
ing values of the modulation amplitude A. In Fig. 5(a), 
the array is straight [A = 0, corresponding to point C of 
Fig. 3(a)] and the wave packet propagates weakly diffract- 
ing with a nonvanishing transverse group velocity along a 
direction which is determined by the normal to the band 



dispersion curve at qa = 7r/4. In Fig. 5(b), the amplitude 
of the modulation is increased to A — 30 /im, correspond- 
ing to &o ~ 4 [see Eq.(6)], i.e. to point A of Fig. 3(a). In 
this case, after the ac field is switched off at z = A ~ 6.7 
mm, the wave packet does not basically break down, how- 
ever it is now refracted in an opposite direction as com- 
pared to the case of Fig. 5(a). This is the clear signa- 
ture that the wave packet is mostly composed by Bloch 
modes of the second band, i.e. that an almost complete 
transition from negative- to positive-energy states has oc- 
curred, according to the predictions of the two- level anal- 
ysis [see point A of Fig. 3(a)]. In Fig. 5(c), the amplitude 
of the modulation is further increased to A = 45 /jm, 
corresponding to $o — 6. According to the prediction 
of the two-level analysis [see point B of Fig. 3(a)], after 
the ac field is switched off the transition probability for 
the wave packet components is ~ 45%. Accordingly, the 
wave packet breaks up at propagation distances z > A, 
with the two nearly balanced wave packets belonging to 
the first and second array minibands refracting along two 
opposite directions. 



IV. CONCLUSIONS 

Vacuum decay and electron-positron pair production 
due to the instability of the QED vacuum in an external 
oscillating field are remarkable predictions of the Dirac 
theory. Their experimental demonstration, however, re- 
quires ultrastrong and high-frequency laser fields which 
are not yet available. In this work a classic wave op- 
tics analogue of vacuum decay has been proposed which 
is based on the formal analogy between the temporal 
dynamics of the one-dimensional Dirac equation for a 
single particle in an external spatially-homogeneous os- 
cillating field and the spatial propagation of monochro- 
matic light waves in a periodically-curved binary waveg- 
uide array. Our photonic analogue enables a simple and 
experimentally-accessible visualization in space of the 
process of pair production as break up of an initially 
negative-energy Gaussian wave packet, representing an 
electron in the Dirac sea, under the influence of an oscil- 
lating electric field. 



Appendix A: Two-level description of pair 
production in the dipole approximation 

In this Appendix we briefly review the two-level de- 
scription of pair production, induced by an oscillating 
field in the dipole approximation, for the one-dimensional 
single-particle Dirac equation given by Eq.(7) in the text. 
In making such an analysis, we introduce the variable 
transformation defined by Eqs.(9), and re- write Eq.(7) 
into the more familiar quantum mechanical form 

dip 

ih— = (cp x - eA x )aip + mc 2 ^, (Al) 
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where p x — —ihd x . The independence of the vectorial po- 
tential A x on the spatial coordinate x corresponds to the 
electric dipole approximation and to the interaction with 
a spatially-homogeneous oscillating electric field. Owing 
to the pure time dependence of the external potential, the 
momentum is conserved during the transition. Therefore 
only transitions between negative- and positive-energy 
states with the same momentum are permitted in the 
pair creation process. Indicating by p — hk the momen- 
tum of the particle at initial time, a solution to Eq.(7) 
with definite momentum is given by 



"4>(x,t) = ip(t) exp(ifcx) 



(A2) 



where the two-component spinor tp = (<p\,fi) T satisfies 
the equation 



(A3) 



The time-dependent elements of the 2x2 matrix Ai in 
Eq.(A3) are given by 



M12 



-M 



22 



M 



21 



kc ■ 



eA x (t) 
h ' 



(A4) 
(A5) 



To investigate field-induced transitions, let us project the 
spinor ip(t) on the basis of the two eigenstates ip± of the 
free electron (i.e., in the absence of the external field) cor- 
responding to the negative and positive energy branches, 
i.e. let us set 



<p(t) = r_(i)<p_ +r+(t)<p+ 



(A6) 



where 



y / p 2 c 2 + [he(p)Tmc 2 } 2 V ±ht (P) ~ mc 



cp 



,(A7) 



he(p) = ^/p 2 c 2 + (mc 2 ) 2 , and r±(t) are the occupation 
amplitudes of negative and positive energy states at time 
t. The amplitudes r_(i) and r + (t) then satisfy the cou- 
pled equations of driven two-level systems 



d 
i— 

dt V r + 



= Z{t) 



r+ 



(A8) 



with matrix elements given by 



Zu = -Z 22 = e(p) 



pceA x (t) 
h 2 e(p) 



Zio — Z 2 l = — 



mc 2 eA x (t) 
H 2 e(p) ' 



(A9) 
(A10) 



Equations (A8-A10) are analogous to Eqs. (17-19) given 
in text for the driven tight-binding lattice (5). In par- 
ticular, according to the approximations used to derive 
the Dirac equation (7) in Sec. II. B, the expressions of 
the matrix coefficients Z ik (t), defined by Eqs. (18) and 
(19), reduce to those given by Eqs.(A9) and (A10), after 
the replacement qa — » n/2 + k/2 [see Eq.(ll)], assuming 
\k\ <C 7r/2, and using the transformation of variables (9) 
that connect the quantum and optical descriptions. 
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